Effects of Coulomb interactions on the splitting of luminescence lines 
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We study the splitting between the right-hand and left-hand circularly polarized luminescence 
lines in a quantum dot under relatively weak confinement regime and resonant high-power excita- 
tion. When the dot is populated with an even number of electron-hole pairs (biexciton and higher 
excitations), the splitting measures basically the Zeeman energy. However, in the odd number of 
pairs case, we have, in addition to the Zeeman and Overhauser shifts, a contribution to the splitting 
coming from Coulomb interactions. This contribution is of the order of a few meV, and shows 
distinct signatures of shell- filling in the quantum dot. 
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I. INTRODUCTION 

This paper is aimed at a theoretical description of the 
splitting of luminescence lines coming from the deexcita- 
tion of a quantum dot in the presence of a magnetic field. 
There are a few factors contributing to this splitting. 

First, there is a contribution coming from the Zeeman 
energy, which is proportional to the magnetic field, B. 
As luminescence lines arise from exciton recombination, 
both the electron and hole Lande factors enter the ex- 
pression for the Zeeman shift. In a nanostructure, Lande 
factors exhibit a dependence not only on the semicon- 
ductor band structure, but also on the geometry. In the 
recent past, for example, extensive studies on the depen- 
dence on well widthi have been conduced. For magnetic 
fields between 1 and 3 Teslas, the typical values we will 
encounter for the Zeeman shifts are 0.05 - 0.2 meV. 

A second contribution to the splitting is the so called 
Overhauser shifts. It is originated in the hyperfine in- 
teraction between the electron and nuclear spins. The 
electron motion around the nanostructure induces a nu- 
clear spin polarization. Once the nuclear polarization is 
set up, the interaction of the electron with the thousands 
of nuclei conforming the dot leads to a measurable shift 
in the luminescence lines. In a first approximation, the 
Overhauser shift does not depend on the magnetic field 
nor on the quantum dot size, but increases up to a sat- 
uration value with the excitation power. The latter fact 
can be understood in terms of the increase of the nuclear 
polarization as the rate of electron-hole pair creation is 
raised. Typical values of the Overhauser shift are around 
0.1 meV. Recent measurements of this contribution were 
done on quantum dots formed from width fluctuations in 
a well 3 -, and on self-assembled quantum dotsi 

There is a third contribution to the splitting, which 
shall appear when the excitation power is high enough 
to create multiexcitonic complexes. Recombination of 
one electron-hole pair from a complex involves Coulomb 
interactions in a nontrivial way. The emission of a right- 
hand circularly polarized photon proves to be not equiv- 
alent to the emission of a left-hand circularly polarized 
one. Typical values for this shift are around 3 meV, i.e., 



the values of energy separations between electronic states 
in the complex. These Coulomb interaction effects have 
not been studied experimentally nor theoretically so far. 
They constitute the focus of our attention. 

The plan of the paper is as follows. In Sec. [n] we 
present a brief resume of k now n ideas explaining the ex- 
perimental results of Refs. |3l4j . which correspond to low 
excitation power and independent exciton recombination. 
This section is intended as an introduction to Sec. IIIII 
where higher excitations are considered. Multiexcitonic 
complexes are treated in the framework of a corrected 
BCS scheme, which takes a proper account of mean field 
Coulomb interactions and Fermi statistics. We compute 
the position and intensities of luminescence lines coming 
from the decay of systems with up to 21 pairs. The main 
results are the absence of Overhauser or Coulomb inter- 
action effects in complexes with an even number of pairs, 
and the shell filling effects in the shifts due to Coulomb 
interactions for systems with an odd number of pairs. 
Concluding remarks are given in the last section. 



II. THE INDEPENDENT-EXCITON 
APPROXIMATION TO THE a+ - a SPLITTING 

The experiments^^ are usually performed in the so 
called Faraday configuration, where the excitation laser 
beam is parallel to the applied magnetic field, and 
backscattering geometry. A schematic view is presented 
in Fig. □ 

The model parameters are choosen to fit (at least quali- 
tatively) the experimental setup described in paper |5j] . A 
quasi twodimcnsional quantum dot is formed from width 
fluctuations in a 4.2 nm-wide GaAs-AlGaAs quantum 
well. The dot area is around 10 3 nm 2 . Taking into ac- 
count that the GaAs lattice constant is ~ .3 nm, the 
number of nuclei in the dot is ~ 10 5 . 

A simplified two-band structure for GaAs is assumed. 
A schematics of the band structure is drawn in Fig. |3 
where electron - hole coupling for each light polariza- 
tion, and Zeeman splitting in each sub-band are indi- 
cated. The effective electron and heavy-hole Lande fac- 
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FIG. 1: Geometry of the experiment. 



tors for the 4.2 nm-wide quantum well are expected to be 
g e « 0.1, gh ~ — 1.2. 1 The Zeeman energy of an electron- 
hole pair is computed from: 



£Ze 



fi B B(g e S e - g h S h ) = HBB(g e + g h )S e , (1) 



where hb — 0.057 meV/Teslas is the atomic Bohr mag- 
neton, and S e , Sh = ±1/2 are the electron and hole spin 
projections over B, which acts along the dot normal. 
g e + gh ks — 1.1 is the exciton Lande factor. 

The quantum dot is resonantly excited, i.e., the laser 
excitation energy is only a few meV above the lumines- 
cence lines. Consequently, in luminescence there is mem- 
ory of how the dot is pumped. Notice that the sense 
of rotation of the electric field is inverted for the emit- 
ted light. This means that by pumping with c+ light, 
for example, we are in fact populating the cr_ mode for 
backward emission. 

As mentioned before, there are around 10 5 nuclei in 
our GaAs dot. Half of them are 69 Ga and 71 Ga nuclei 
with relative abundances 60% and 40% respectively. The 
second half corresponds to 75 As nuclei. For all these nu- 
clei, the total spin is I = 3/2. Their magnetic moments 
(in units of the nuclear magneton, fi n ) are, respectively, 
2.017, 2.562 and 1.439A 

The hyperfinc interaction between one conduction 
band electron and nuclei in the quantum dot is described 
by the Hamiltonian: 



"Hen = a^S e - I n \ 1p e (f n ) | 2 



(2) 



where a = (2/3)/io5effn/^sMn, and /io is the magnetic 
permitivity of vacuum. ip e is the electron wave function, 
which contains both an envelope (orbital) part and the 
Bloch functions, which are s-functions. tp e is evaluated 
at the positions of the nuclei, i.e. the origin of each cell. 
For Bloch s-functions, 0, we have |</>(0)| 2 ~ Z 3 / (a^n) 3 , 
where Z is the atomic number (31 for Ga and 33 for As), 
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FIG. 2: Schematics of the assumed band sructure in the dot. 
The electron-hole pairs created in the absorption of <r+ and 
(j— photons are drawn. 



as is the Bohr radius, and n = 4 corresponds to the 
4s electrons conforming the conduction band. The fact 
that the orbital wave function is extended all over the dot 
(meaning that the square of the orbital wave function is 
proportional to 1/N n , where N n is the number of nuclei 
in the dot) make the hyperfine contribution to the energy, 
Eq. (J2J), roughly independent of the dot size. In a mean 
field approximation, (J) is oriented along B, and we get 
for the interaction energy of one electron with the 10 5 
nuclei in the dot: 



^hyperfine — Sen (I)Se 



(3) 



The g en constant for bulk GaAs was found to be around - 
0.090 rneVi A similar value was observed in the quantum 
dots studied in Refs. and|jj 

Finally, there is a contribution to the exciton energy 
coming from the electron-hole exchange interaction in 
a deformed quantum dot. For the quantum dot under 
study 8 , we have 



(exch) 
-± 



=F 0.012meF. 



(4) 



With the help of Eqs. llll.'il II) . and the independent 
exciton approximation, we write the single pair energy 
in the form: 
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e±=e x + £± + VBB(g e + g h )S e + g en (I)S e , (5) 

where the + index corresponds to S e = 1/2 (er+ polariza- 
tion of the emitted light), and the - index to S e = —1/2 
(<r_ polarization of the emitted light), ex is the contri- 
bution of the electron-hole Coulomb attraction. e± gives, 
respectively, the position of the tr + and cr_ luminescence 
lines. The shift between them is 

Ae = £_-£+ = e (exc ^ - ^ B B(g e + g h ) - g en (I). (6) 

It depends on the excitonic populations, N+, N-, 
through (I). The Overhauser shift is, by definition, 

^^Overhauser — £/en (-0 ■ (7) 

It can be measured by turning off the contribution of each 
nuclear species to the mean nuclear spin with the help of 
microwave radiation sintonized to the frequency of the 
corresponding nuclear magnetic resonance^. Under a + 
pumping, the electron spin polarization is (S e ) < 0, and 
the mean nuclear spin is (/) > in order to minimize 
the hyperfine interaction energy. This means that the 
Overhauser shift is added to the Zeeman splitting. Under 
<r_ pumping, on the contrary, the Overhauser shift is 
substracted from the Zeeman splitting. 



III. ACCOUNT OF MANY-PARTICLE 
COULOMB INTERACTIONS 

In the independent-exciton approximation, Eq. JSJ, 
the dependence on pumping comes only from (/). In 
the experiments reported in Refs. I'l ll the mean exciton 
number in the dot is less than one. In paper 0, there is a 
study of the dependence of Ae on the excitation power, 
but what is indeed observed in the single-exciton line 
is the saturation of the nuclear polarization, (/), as the 
occupation of the exciton state increases. 

In the present section, we shall discuss the more general 
situation, where the excitation power is enough to create 
a few pairs in the dot. Many-body Coulomb interactions 
may introduce an additional dependence of Ae on iV + 
and N- . The energy of a system with N+ , N- pairs is 
written as: 



E(N+,N_) = E C oui{N+.N_) 

+ v B B(g e +g h )(N + -N_)/2 

+ g en (I)(N+-N-)/2, (8) 

where EqouI accounts for Coulomb interactions among 
particles. The positions of the er + and <j_ lines are ob- 
tained from: 



s+ = E Co ui{N + ,NJ)-E C oui{N + -l,N_) 

+ [i B B(g e +g h )/2 + g en {I)/2, (9) 



e_ = E C oui{N +: N-) - E C oui(N+, iV_ - 1) 

- fi B B(g e +g h )/2-g en (I)/2. (10) 

Coulomb interactions introduce a contribution to the 
shift: 



Aecoul = E Co ul(N + -l,N_)-Ecoul(N + ,N--l). (11) 

It is reasonable to accept that Ecoui does not vary 
when N + and iV_ are interchanged: Ecoui(N+, N-) = 
Ecoui{N-,N+). 

In order to compute Ecoui , we use a model of parabolic 
dot with harmonic confinement, fku w 10 meV. This value 
of to is chosen to reproduce the exciton diamagnetic shift 
of 0.025 meV/Teslas 2 reported in Refs. 00. The corre- 
sponding oscillator lenght is around 12 nm, a value simi- 
lar to the dot characteristic dimensions. Effective masses 
m e — 0.067 mo, = 0.15 mo, and relative dielectric 
constant k — 12.5 are used. Electron- hole pairing is ac- 
counted for by means of a BCS wave function. Notice 
that the Bohr radius for GaAs is aroud 7 nm. This means 
that the area occupied by an exciton is around 150 nm 2 . 
In our 10 3 nm 2 dot, 7 excitons are already closed packed, 
and the effects of Fermi statistics should be important. 
These effects and the mean field Coulomb interactions 
are correctly described by the BCS function. 

Fluctuations in the particle number, not conserved 
by the BCS function, are important for a system with 
around 10 pairs. Thus, the BCS function should be pro- 
jected onto a subspace with fixed particle number. The 
projection scheme we use is the so called Lipkin-Nogami 
schem o 10 ' 11 ' . In the present situation, we shall use a 
Lipkin-Nogami scheme with two conserved charges, N + 
and N-. As this situation is not common in the litera- 
ture, we give in the Appendix a brief description of the 
method and the explicit equations to be used. 

Let us first discuss the simplest case beyond the inde- 
pendent exciton approximation: the cr+ — ct_ splitting in 
the biexciton line. 

The lowest biexciton state has (N+, N-) = (1,1). The 
first state with N + or 7V_ equal to 2 has excitation en- 
ergy higher than 8 meV (with our model parameters). 
Under quasiresonant excitation and temperatures below 
4 K, only the ground state will be populated with high 
probability. 

The Coulomb contribution to the splitting, Eq. ^ L 1 p . is 
thus zero for the biexciton lines. It is not difficult to see 
that the Overhauser shift is also zero or, at least, much 
smaller than the shift for the exciton. Indeed, an intuitive 
argument suggests that, as the spin polarization is zero in 
the ground state, the mean nuclear spin will be zero too. 
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FIG. 3: (Color online) (a) Coulomb contribution to the <r+ — 
a- splitting in the (N+,N-) = (2, 1) system, (b) Intensities 
of the luminescence lines at B = 1 and 3 Teslas. The reference 
energy is E gap - Only Coulomb interactions are included when 
computing the line positions. 



In fact, the situation is a bit more complicated. One can 
guess that, as each excitation has a certain weight in the 
density matrix, exciton states, including the long-lived 
dark states, play an important role in determining the 
mean nuclear spin 3 . Nevertheless, a lower value for (I) 
is expected and, consequently, a lower shift due to the 
Overhauser effect. 

There is, nevertheless, a very important effect of 
Coulomb interactions consisting in the concentration of 
the luminescence in a single, coherent, line, as will be- 
come clear below. 

The conclusion is that the biexciton lines exhibit a 
weaken Overhauser effect, and no shift coming from 
Coulomb interactions among particles. This conclusion 
can be extended to any states with N + = N-, i.e., when 
there is an even number of pairs in the dot. 

The situation is different, however, when the num- 
ber of pairs is odd. Let us consider, for example, the 
(N + ,N-) = (2,1) case. We expect a nonzero mean 
nuclear polarization and, consequently, a nonzero Over- 
hauser shift. The magnitude of this shift is expected to 
be around 0.1 meV. The Coulomb contribution to the 
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FIG. 4: (Color online) (a) Splitting induced by Coulomb in- 
teractions as a function of the number of pairs in the dot, 
(b) Comparison between single-pair and coherent cr+ lumi- 
nescence in the (N+,N-) = (7,6) system at B = 1 Tesla. 



splitting is much greater, of the order of a few meV. 

We show in Fig. HJa) the Coulomb splitting for the 
luminescence lines coming from the decay of the (2,1) 
system as a function of the magnetic field, B. Notice 
the smooth dependence of Aecoui on B. We shall stress 
that, when computing the position of the lines, e±, we 
assume transitions from the ground state of the ( N + , iV_ ) 
system to the ground states of the (N+ — l,iV_) and 
(N+,N- — 1) systems. We verified that these so called 
coherent luminescence lines account for more than 80 % 
of the total luminescence (see Appendix). We represent 
in Fig. EJb) the relative intensities of both lines at B = 1 
and 3 Teslas. Only Coulomb interactions are included 
in the line positions. The reference energy is the band 
gap. Lorentzians with a width T = 0.2 meV are used to 
represent the lines. The variation of the magnetic field 
induces a blucshift on both line positions, which is of 0.45 
meV for the a + line, and 0.22 meV for the er_ line. For 
comparison, let us notice that the Zeeman energy (not 
included in the figure) moves the <t+ line 0.062 meV to 
the left, and the <r_ line 0.062 meV to the right, when B 
varies from 1 to 3 Teslas. 

The luminescence lines coming from the decay of the 
(2,1) system shall be observed under er_ excitation. Due 
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to the reasons mentioned above for the biexciton, the 
excitation of the (3,0) system is hardly possible. Under 
<7_|_ excitation, Aecoui changes sign, and the lines inter- 
change the intensities. 

It is interesting to look at still higher excitations. To 
be definite, we consider cr_ pumping and follow only the 
luminescence lines coming from the decay of the (n+1, n) 
systems. The splitting induced by Coulomb interactions 
is drawn in Fig. EJa) as a function of the total number 
of pairs, N = 2n + 1. The splitting is close to the har- 
monic confinement energy, Hu>o = 10 mcV, for the N = 3 
system, but only around 3 meV for the N=21 complex. 
We observe also distinct signatures of shell filling at the 
magic numbers N — 1 = 2, 6, 12 and 20. 

The effects of Coulomb interactions are made evident 
in Fig^Jb), where single-pair and coherent cr + lumines- 
cence in the (N + ,N_) — (7,6) system at B = 1 Tesla 
are compared. The single-pair luminescence is evaluated 
from the pair energies and occupations. We see the two 
main effects of Coulomb interactions in this figure. First, 
there is a strong enhancement of luminescence intensity. 
This is the analog of the well known exciton effect. Sec- 
ond, the coherent peak is shifted with respect to the po- 
sition of the single-pair peaks. 

Finally, we shall stress that, for a multiexciton system 
with TV > 2 to live in a quantum dot, the confining po- 
tential should correspond to a deep well, because these 
systems are, most likely, unbound^. This does not seem 
to be the situation reported in [50, where the dot is 
formed as a result of fluctuations in the well width. It 
does not represent, however, a limitation to the validity 
of our results, because quantum dots with the desired 
properties can presently be grown at will. 



IV. CONCLUDING REMARKS 

We have shown that luminescence lines arising from 
the decay of one electron-hole pair in a multiexcitonic 
complex may exhibit a strong asymmetry with respect 
to the sense of circular polarization of the emitted light 
if the number of pairs in the complex is odd. Let us stress 
the conditions under which this statement is expected to 
be valid. 

First, it was assumed that N + — N- = ±1. Lumines- 
cence lines from the complex with N + = N- were shown 
to exhibit roughly no Overhauser nor Coulomb shifts. 
From an experimental point of view, a small disbalance 
between N + and iV_ is expected under quasiresonant ex- 
citation and very low temperatures. High excitation pow- 
ers are needed to create the multiexcitonic complexes. 

A second, very important, assumption is the coherent 
(collective) character of the decay process. Lines which 
can be interpreted as single- or independent-pair decays 
are possible too. One can guess that these single-pair de- 
cays should play an important role in luminescence un- 
der a strong confinement regime and well above band gap 
excitation^. In our calculations, we assumed an inter- 



mediate confinement regime, Hujq — 10 meV, in which 
case effects due to Coulomb interactions become very 
important, and the coherent luminescence accounts for 
more than 80 % of the total luminescence. The effects of 
Coulomb interactions were shown to be basically the fol- 
lowing: and enhancement of luminescence intensity and 
concentration on a single line, which is shifted with re- 
spect to the position of the single-pair decays. Band- 
filling in the quantum dot is seen as a blueshift of the 
peak position. 
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APPENDIX A: LIPKIN-NOGAMI SCHEME 
WITH TWO CONSERVED CHARGES 

The Hamiltonian describing the electron-hole system 
in the two-dimensional quantum dot is written as: 

H = E( < « )e « e «+ < « i)/i « /i «) 

n 

+ g ^2( r ,s\l/r\u,v) e\e\e v e u 

rsuv 

+ — ^ (r, s\l/r\u,v) hlh\h v h u 

rsuv 

- p"^2(r,s\l/r\u,v) elh%e u , (Al) 

rsuv 

where, for simplicity, we have not included electron-hole 
exchange terms. j3 = 0.8 e 2 / (Aireonlo) is the strength of 
electron-electron Coulomb interactions, whereas j3' and 
j3" characterize the hole-hole and electron-hole interac- 
tions. The calculations presented in the paper correspond 
to the symmetric situation, in which j3 — (3' — f3" . A fac- 
tor of 0.8 is included in (3 to approximately account for 
the quasi two-dimensionality of the dot. Iq = y/h/(m e £i) 
is the oscillator length corresponding to the modified fre- 
quency £1 (to be discussed below). The single-particle 
electron and hole energies in a harmonic oscillator po- 
tential and a magnetic field are given by: 

4 e) =Ml (2k + \l\ + l) + fkj c l/2, (A2) 



t { n h) = — {tM (2k + |/| + 1) - hu c l/2} , (A3) 
mh 

where = \Joj 2 -j-tjf/4, and uj c = eB/m e is the electron 
cyclotron frequency. It is assumed that the confinement 
potential for holes is such that its characteristic length is 
l also. 
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n,r,s, etc are two-dimensional harmonic oscillator 
states with characteristic length l . They are specified 
by the radial quantum number, k, the orbital quantum 
number (angular momentum projection over B), I, and 
the spin projection, a. Notice that we have not included 
the Zeeman and Ovcrhauser contributions to the single- 
particle energies (JA2IA3(1 . Notice also that the Coulomb 
matriz element, (r, s|l/r|w, v) is zero unless the spin pro- 
jections satisfy: a r — a u , a s — a v . 

There are four conserved charges commuting with the 
Hamiltonian, (|A1|I . They are: 



N. 



e1 



A^ 



el 



(A4) 



N, 



hi 



h i h * 



N, 



I: I 



where the notation means that only spin-up electron 
states enter the sum in Af e j, etc. 

In the BCS function we shall write, the number of spin- 
up electrons exactly coincide with the number of spin- 
down holes, and the number of spin-down electrons with 
the number of spin-up holes. Let us define A^ + = N e -\, 
and Af_ = iV e j_. The mean value of the Hamiltonian, 
Eq. IjAlJl . should be minimized with the constraint that 
the mean values (A/+) and (A^_) should be fixed. Mean 
values are obtained by averaging out with the function: 



\BCS) - l[(u r + i; re tftt) Y[(u s + v s el4) |0), (A6) 



where |0) denotes the vacuum state, and the "conjugate" 
hole state n is defined in terms of the electron state 
n = (k,l,a) as n — (k,—l,—a). The coefficients u r ,v r , 
satisfying the normalization conditions u 2 , + v 2 = 1 , are 
to be used as variational parameters. Notice that 



(A7) 



nl 



The standard BCS equations (gap equations) are ob- 
tained by introducing Lagrange multipliers which 
are fixed from the equations i|A7(l . and minimizing the 
function: 



(/3 + f3')/2j2(n, m\l/r\m, n)v 2 n v 2 m 
l3"^(n,m\l/r\m,n)v n u n v m u m ). (A9) 



Notice the exact cancellation of direct Coulomb interac- 
tions in the symmetric case, j3 — j3' = f3" . 

A rather nontrivial fact, which helps understanding the 
Lipkin-Nogami scheme, is that the BCS gap equations 
can also be interpreted as a search for the optimal linear 
in N+ and N- approximation to H . Indeed, let us define 
P = Ao + A + iV + + A_A^_, and require the minimization 
of the mean square deviation: 



(A5) 

The equation: 



G(X) = ((H-P) 2 



dG/dXa = (H-P) = 0, 



(A10) 



(All) 



is used to fix the Ao parameter, whereas the equations 
dG/dX+ = dG/dX- = prove to be equivalent to the 
BCS gap equations. They can be written in the form: 



= ((H-P)N+), 



= ((H - P)N. 



(A12) 



(A13) 



The chemical potentials, fj,±, can be identified with the 
A± parameters in the present case. 

In the Lipkin-Nogami method, we make a step forward 
and look for a quadratic approximation to H: 



P = Ao + A + iV + + X-N- + X++N 2 
+ X+-N+N- + A__A^. 



(A14) 



Minimization of G in Eq. I|A10|) leads, in addition to 
l|AlllA12IA13j) . to the equations: 



0=((H-P)Nl), 



(A15) 



F(v) = (H-p+N+-n-N-), (A8) 



= ((H-P)Ni), (A16) 



with respect to the variational parameters v n . The mean 

value of the Hamiltonian is given by: q _ i/jj _ p\jy ]\j \ (A17) 

V ( (e) ( h A 2 Ec l s - HA11IA12IA13IA15IA16IA17|) conform a linear 

' ' — 2-^ ~^tfi J v n system from which we obtain the A parameters in terms 

n of averages like (HN+), (N+), etc. In fact, we need only 

+ (/? + /?' — 2f3")/2 / (n, m\l/r\n, m)v 2 v 2 ^ explicit expressions for A ++ , A H and A : 
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A+_ = 



(HN+N_) - {HN+)(N_) - (N+)(HN_) + (H)(N+)(N_) 



(Nl)(Nl) - (JV|>(2V_> 2 - <VV+> 2 (A^) + (iV + )2(7V_)2 > 
((gjV 2 ) - (H)(Nl))((Nl) (iV + ) 2 ) - - (H)(N + ))({Nl) (A + )(7V 2 )) 

I 



(A18) 
(A19) 



and similarly for A The parameters A+ and A_ are 

absorbed into the definition of the chemical potentials: 



/ i + = A + + 2A ++ <7V + ) + A + _(iV-}, 



(A20) 



(and similarly for fj,—). They are obtained from Eqs. 
(|A7|) . in which we introduce the standard parametriza- 
tion: 



<T = - I 1 - 



£ n1 



(A21) 



and similarly for v 2 ^. The pair Hartree-Fock energy is 
given by: 



HF = .(e) , Ah) _ R n 



n 1 tl 



/3"(n,n|l/r|n,n)-A++«-<) 



+ (/3 + /3' - 2/3") (n,k\l/r\n,k)vl 



[fi + /3' - 2/3") J2(n, k\l/r\n, k)i 



A' I 



- (/3 + /3') ^ (n,k\l/r\k,n)i 



(A22) 



The gap parameters, A+, introduced in Eq. (|A21J| . are 
obtained from Eq. I|A12|I . The latter can be put in the 
form of gap equations: 



A+ = P" s\l/r\s, n) Af 

si J(er- M+ )2 + (A+)2 



(A23) 

and similarly for A~. 

The iterative procedure designed to solve these equa- 
tions is as follows. We start from a set {s^ 1 , A n }. Eqs. 
(|A21|I are introduced into Eqs. I|A7() and the chemical 
potentials are found. Then, we compute A ++ , 

A^ and A Hartree-Fock energies and gap functions 

are recalculated from Eqs. (jA22IA23j) . and the process is 
repeated until convergence is reached. 

For completeness, let us write the explicit expressions 

for the mean values needed to compute A++, A^ and 

A__ : 



= (N + y + (N + )-J2vt 



(A24) 



(N*) = (N + ) 3 + 3(N + ) 2 + (N + ) 

- 3 ((N + ) + l)$>n + 2£«n, (A25) 



n1 



"T 



(Nt) = (7V+) 4 + 6(AT+) 3 + 7(7V+) 2 + {N + - (6(N + f + 18{N + ) + 7) ]T v* 

+ 4(2(7V + ) + 3)^^ + 3(^ V 4 J 2 -6^^, (A26) 

nf nf nt 

i 

and similarly for (iV 2 ), etc. Notice that, because of the factorize, (N+N-) = (N+) (N-) . Concerning mean val- 
factorizable form Eq. !A6jl . mean values like (N+N-) ues of H with powers of N+, we have: 



(HN+) = (H)(N + ) + ^(4 e) +i h) )v 2 n u 2 n 
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+ {(P + V -W")(n,m\l/r\n,m)-((3 + (3')(n, m\l/r\m,n)}v 2 n v 2 m (l - v 2 n ) 

- P" X! (n,m\l/r\m,n) v n u n v m u m {\ - 2v 2 ), (A27) 

nf,mf 

(HN+N-) = (HN+)(N-) + (HN^)(N + ) - (H)(N+)(N-) 

+ (/? + /?'- 2/3") ]T (n,m\l/r\n,m) v 2 n u 2 n v 2 m u 2 ml (A28) 

(HN 2 ) = {H){N 2 )+2{N + )((HN + )-(H)(N + )) + Y,(ti ) +t { n ) )u 2 l vl(l-2v 2 n ) 

"T 

+ (P + P' -2f3")J2(^^/r\n,m)v 2 n v 2 m u 2 n (u 2 n -v 2 l ) 

nf ,m 

+ (/3 + ~ 2/3") X! (n,m\l/r\n,m) v 2 n v 2 m u 2 n u 2 m 

n|,mf 

- + ^ (n,m\l/r\m,n)(2v 2 n v 2 m -5vy m + 2vy m + vtvt) 
- /?" ^ {n 1 m\l/r\m,n)v n u n v m u m {l-&v 2 l + Av i n + 2v 2 n v 2 m ). (A29) 

i 



Notice that, from Eq. (|A28|) it follows that A+_ = 
in the symmetric system with (3 = (3' = j3" . Finally, 
the Lipkin-Nogami energy is obtained by replacing the 
operators N + , A_ in P by its eigenvalues. It can be 
written in the form: 

E LN = (H) - \ ++ ({N 2 + ) - {N + ) 2 ) 

- X—((N 2 ) - (N-) 2 ). (A30) 

To get an idea of the intensities of the luminescence 



lines, we compute the matrix elements of the interband 
dipole operator: 

D+=Y,enhn. (A31) 

"T 

A similar expression can be written for D_ . The intensity 
of the transition from the ground state of the ( N + , A_ ) 
system to the ground state of the (N + — 1, iV_) system 
is proportional to the matrix element of D + squared: 



I+ = \(N + -l,N-\D + \N+,N-)f 



Y[(u' n U n +v' n V n ) 



•>.[ 



^2u' n v n Y[ {u' n u n + v' n v n ) 



(A32) 



where (u,v) and (u',v') are, respectively, the BCS coef- 
ficients for the initial and final states. It can be verified 
that the transition from the ground state of the (N+, iV_) 
system to the ground state of the (N + — l,iV_) system 
accounts for most of the transition probability. Indeed, 
a total intensity can be defined in the following way: 



jtotai = YMD + \N + ,N^)\ 2 



= E^ + (£«»«*) 2 , (A33) 

where the sum runs over all possible final states, 
For the (2,1) complex at B = 1 Tesla, with the parame- 
ters used in this paper, I+/I+ tal = 0.84, whereas for the 
(10,9) system, I+/I^ tal = 0.92 . 

The two terms in Eq. I|A33(1 have a simple interpreta- 
tion. If pairing is neglected, the decay of a pair in the 
state n proceeds independently of the rest of the states. 
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The matrix element of D + is equal to the probability 
of finding a pair in the state n, i.e., to v\. Thus, D + 
squared equals v\. The first term in Eq. (|A33j) is hence 
the probability of decay of individual pairs. The sec- 



ond term comes from pairing. In our model, at B = 1 
Tesla, I i ^' tal from the (7,6) complex, for example, has a 
component coming from individual decays equal to 5.20, 
whereas the pairing contribution reaches the value 59.34. 
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